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THE DYNAMICS OF AN ARTICULATED n-TRAILER VEHICLE 


ALEJANDRO BRAVO-DODDOLI AND LUIS C. GARCIA-NARANJO 


Abstract. We derive the reduced equations of motion for an articulated n-trailer vehicle that moves under its own 
inertia on the plane. We show that the energy level surfaces in the reduced space are (n + l)-tori and we classify the 
equilibria within them, determining their stability. A thorough description of the dynamics is given in the case n = 1. 


1. Introduction 


We consider the dynamics of an articnlated n-trailer vehicle that moves nnder its own inertia. 
Snch system consists of a leading car, or trnck, that is pnlling n trailers, like a Inggage carrier in the 
airport. The leading car and the trailers form a convoy that is snbjected to (n -|- l)-nonholonomic 
constraints, one for each body. 


This system is a canonical example in nonholonomic motion planning, which is fnndamental in 
robotics, and has been extensively considered from the control perspective (see e.g. [I21EIII3] and 
the references therein). 

The constraint distribntion defined by the n-trailer system has also received interest in differential 
geometry. It is a Gonrsat distribntion (see 0) and, as shown in na, all possible Gonrsat germs of 
corank n -|- 1 are realized at its different points. 

The dynamics of wheeled vehicles moving on the plane has been considered in e.g. [h El |3]. The 
hrst two of these deal with certain properties of the n-trailer system. However, the majority of the 
existing references to this system deal with its kinematics and disregard its dynamical aspects. To 
the best of onr knowledge, a detailed stndy of the dynamics of the n-trailer system is missing. 


The nonholonomic constraints in the n-trailer system arise by assnming that each of the bodies in 
the convoy has a pair of wheels that prohibit motion in the direction perpendicular to them. Each 
of these constraints is identical to the nonholonomic constraint for the well-known Ghaplygin sleigh 
problem [5]. Hence, the articulated n-trailer vehicle is a generalization of the Ghaplygin sleigh system 
that is recovered when the number of trailers n = 0. (Other generalizations of the Ghaplygin sleigh 
are considered in 0 )- 


We mention the thorough understanding of the dynamics of the Ghaplygin sleigh has resulted in 
the design of control algorithms, where the control mechanism moves the center of mass [T7]. We 
are hopeful that the results in this paper will turn out to be useful for control purposes. 


The outline of the paper is as follows. In Section we introduce the system and the notation 
that we will follow in our work. We define the configuration space, the nonholonomic constraints 
and we write down the kinetic energy of the system. We also discuss the related SE{2) symmetry 
of the problem and describe the reduced space. In Section we derive the reduced equations of 
motion (3.1). Our approach follows the method suggested in [8]. We show that the energy level 
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sets are diffeomorphic to {n + l)-tori, and we give working expressions for the restriction of the flow 
to them. Section considers the eqnilibria of the system assnming that the center of mass of the 
leading car is displaced a distance a > 0 from the wheel’s axis. We give a complete classihcation of 
all the eqnilibria in an energy level set and perform their linear stability analysis. It is fonnd that the 
straight line motion of the convoy in the direction of the center of mass of the leading car and with 
all of the trailers aligned behind it, is asymptotically stable. In Section]^ we deal with the case where 
the center of mass of the car lies on its wheels’ axis. We give necessary and snfficient conditions on 
the velocities for the existence of eqnilibria of the rednced system and we do an exhanstive treatment 
of the case n = 1. Finally, in Section we comment on the interest to analyze the inflnence of the 
singnlar conhgnrations on the dynamics. 


2. The n-xRAiLER mobile vehicle 


Following and other references given in these works, we consider a mnlti-body car system 

(Bq, Si, ..., 'Bn) that consists of a car Bq pnlling n trailers. Si,..., S„. The trailers form a convoy 
(like in a Inggage carrier) that moves on the plane (see Fignre 2.1 for the case n = 2). 

Each body in the convoy has a set of wheels and we denote by {xi,yi) the coordinates of the 
midpoint of the wheel’s axis {i = 0,... ,n) with respect to a chosen cartesian frame. The orientation 
of Sj is determined by the angle 6i between the main axis of the body and the x axis of the chosen 
frame (see Fignre 2.1). 


2.1. Kinematics. The convoy condition reqnires that the body Sj is hooked to the preceding body 
Sj_i. Following [T^ [13] we assnme that the hook ing is done via a link of length £ that connects 
{xi,yi) with (xj_i, as illnstrated in Fignre 2.l|p The hooking of the convoy thns dehnes the 2n 
holonomic constraints 


Xi+ £ cos 9i — Xi_i = 0, yi + £sm9i — yi_i = 0, i = (2.1) 


On the other hand, the wheels on each of the cars impose a nonholonomic constraint that forbids 
any motion of the given body in the direction perpendicular to its main axis. In this way we get the 
n + 1 nonholonomic constraints 


Xi sin 9i — iji cos = 0, z = 0,1, 


.n. 


( 2 , 2 ) 


In view of the holonomic constraints (2.1), the conhguration of the convoy is fully determined by 
the value of the coordinates 

^ 0 , Vo, 9o,9i,..., 9n- 

Therefore, the conhguration space of the system is the n + 3 dimensional manifold Q = SE{2) x T” 
where SE{2) denotes the Euclidean group in the plane and T" is the n-torus. The nonholonomic 
constraints (2.2) dehne a rank 2 constraint distribution D on Q. 


^Other hooking mechanisms are possible and have been considered in the literature, 
various models, including the one that we consider in this paper m- 
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Figure 2.1. The n-trailer vehicle with n = 2. 


2.2. Dynamics. We assume that the center of mass of the leading car Sq is displaced a distance a 


from the midpoint of its wheel’s axis (xo,t/o) along the principal axis of the body (see Figure 2.1). 
Therefore, if {xc,yc) denote the coordinates of the center of mass of 'Bq, we have 

xc = xo + acos9o, yc = yo + asm9o. (2.3) 

We will denote the total mass of Bq by M and its moment of inertia about its center of mass by 
Jq. On the other hand, we shall suppose that the trailers “Bi,, Bn are identical, with their center 
of mass lying at the midpoint of the wheel’s axis {xi,yi). Their total mass is denoted by m and the 
moment of inertia about {xi,yi) by J. 

The kinetic energy of the system is given by 


V = 2 I ^o95 + + yt;) + 




2 = 1 


n 

2=1 


Using (2.3) we get 


= - j ( Jo + Ma^)9l + M{xl + yl) + 2Ma9o{yo cos 9o - xq sin 6^0) + J ^ + m ^(J- + yi) 

1=1 


2 = 1 
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The Lagrangian of the system L : TQ —)• I 
of the coordinates (a:o, Vo, ■ ■ ■, On) of Q. 
constraints (2.1) imply 


“ is obtained by expressing the above qnantity in terms 
In order to eliminate (xj, i/i) we note that the holonomic 


Xi = Xq- i'^cosOj, yi = yo- i'^sinOj, i = (2.4) 

i=i i=i 

Differentiating the above and adding yields, 

n n n 

+ yl) = n{xl + ijl) + 2£ + 1 - j)9j{yo cos % - xq sin Oj) + f ^(n + 1 - j)e] 

i=l j=l j=l 

n—1 n 

+ 2 ^^ 1 “ 3)0k0j cos{0k - Oj), 

k=l j=k-\-l 


where we have used the identity 


E 

i=l 


^T,cos0n + 


0=1 


0=1 


n n 


'^{n + l-j)Tf + Y^ 2{n + l-j)TkTjCos{ek-ej), 

j=l k=l j=k-\-l 


that holds for arbitrary scalars Ti,..., Tj0 

Therefore, the Lagrangian of the system £ : TQ —)■ M is given by 


- [{Jo + Ma^)0'^ + (M + nm) (x^ + y"^) + 2Ma6{y cos 0 — x sin 6) 

n n 

+2mi + 1 — j)Oj{y cos Oj — x sin Oj) + + (n + 1 — j)mi‘^)0 


(2.5) 


+2mf Y (^ + 1 - j)OkOj cos{Ok - Oj) J , 
k=l j=k-\-l J 

where we have introduced the simplihed notation x = xo,y = yo,0 = Oq. 

Using again ( 2.4[ ), we can write the nonholonomic constraints (2.2) in terms of the coordinates 
{x,y,0,0i, ...,On)oiQas 


X sin Oi — y cos Oi + i'Y^cos{Oi — Oj)Oj = 0, i = 0,...,n. (2.6) 

i=i 


In principle, using (2.5) and (2.6), one could write down the equations of motion for the system 
in terms of Lagrange multipliers using the Lagrange-D’Alembert principle (see e.g. m)- However 
this approach does not make use of the symmetry of the problem that we discuss next. 


2 


We use the convention that a sum over an empty range of indices equals 0. 
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2.3. Symmetries. The system possesses an SE{2) symmetry associated to the arbitrariness of the 
origin and orientation of the chosen cartesian frame. The action of the matrix 


cos (p — sm (p r 

g= I sin(p cos(p s 1 ^ SE{2) 

0 0 1 

on the conhguration {x, y,9,9i,..., 9n) G Q is given by 

g ■ (x, y, 9,9i,..., 9n) = (aio cos 93 — j/o sin(y 9 + r,x sin ip + y cos ip + s,9 + ip,9i + ip,... ,9^ + (p). 


It is immediate to check that the Lagrangian (2.5) and the constraints (2.6) are invariant under 


the tangent lift of this action. It follows that the equations of motion drop to the quotient D/SEi2) 
which is a rank two vector bundle over the n-torus ¥"■. 

We denote the angles between subsequent bodies in the convoy by 


ai = 9 - 9i, ai = 9i_i - 9i, i = 2, 


(2.7) 


see Figure 2.1 The value of these angles is invariant under the SE{2) action dehned above and their 
values serve as coordinates on the base of the reduced space T)/SE{2). 

Next, we denote by u the component of the linear velocity of the leading body 'Bq along its main 
axis, and by oj its angular velocity. We have 


u = X cos9 + ysin9, u = 9. 


As it shall become clear below, the variables u,uj serve as linear coordinates on the hbers of the 
reduced space T)/SE{2). The reduced equations of motion form a set of n + 2 nonlinear, coupled, 
hrst order ordinary differential equations for u, u, ai,..., 


3. The equations of motion 
The purpose of this section is to show the following. 

Theorem 3.1. The reduced equations of motion of the n-trailer vehicle are given by 

1 


u = 

2R(a) 

to = 

Mauu 

Jo + Ma2’ 


Msinoi 

Ol.\ = 

£ ’ 



dfc = 

-j 1 n 


dR 

dai 


2 , < 5 («) , Ma 2 

t^R{a) R{a) 


(3.1) 


0=1 
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where the coefficients Ak are defined by (|3.5|) below and 

j 


Q{a) := cosaisinai j mi'^ E n 

V j=l \k=2 

11 j 


R{a) ■= M + m En 

\j=i k=i 

where we denote a = {ai ,..., 


cos ak 


J 




cos ak 


k=2 


(3.2) 


cos ttfc I + ^ 


+ 7^ ’ 


k=l 


The proof of this theorem follows the approach developed in [S] to obtain the equations of motion 
of regular mechanicaQ nonholonomic system. 


We begin by noting that the relations (2.7) imply 


9i = 9 — aj, i = 1,..., n. 


(3.3) 


Using these expressions, we can write the nonholonomic constraints (2.6) as 
i:sin6* — ycos9 = 0, 


xsm 


^ ~ "i ) - 2/ cos 


9 - ^ aj \ + £ ^ cos I ^ ttfe I I 0 ^ dz I = 0, i = 1,... ,n. 




j=i \k=j+i 


1=1 


nd nd ^ , d 

Zi = cos6*—+ sm0—+ > Ak 

Ht nil < ^ 


Use (x, I/, 0, , arfi) as coordinates on Q and consider the vector fields on Q 

d d 

> 4,. - y.r, = 

dy 

where 

1 \ 

= J COS ajj (sin ak-i - cos ak-i sin a^), A; = 1,..., n. 


(3.4) 


(3.5) 


In the above expression and in the sequel, we use the convention that the product over an empty 
range of indices equals 1 and ao = 0. 

It is readily seen that Zi and Z 2 are linearly independent. Moreover, using the identities 

U-i 


sm 


E 


aj 1 


0=1 
j / 1-2 


E“d E “d n 

j=l \k=j+l / \s=l 


COSOs I sinoj. 


U-i 


( JJ^cosOs I (sina;_i — coso^-i sino;) = — j J^cosOs ) sino^. 


(3.6) 

(3.7) 


1=1 \s=l 




^Note that R(a) > 0 for any value of a. 

“^By regular mechanical we mean that the Lagrangian is the kinetic energy minus the potential energy, where the kinetic energy defines 
a Riemannian metric on the configuration manifold, and the constraint distribution has constant rank. 
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one can check that Z\ belongs to V. It is easy to see that Z 2 is also a section of D. It follows that 
{Zi, Z 2 } is a basis of sections of T> and any tangent vector v E TQ belonging to D can be written as 
a linear combination 


V = uZi + 0 JZ 2 , 

The components of the above equation give 
X = ucos9, 
y = usm9, 

9 = u, 

u sin 


u,u E 


(3.8) 


ai = u 


u 

= Y 


(3.9) 


/'k-2 


i 


COS aj I (sin ak-i — cos ak-i sin Uk) , k = 2,... ,n. 


\j=i 


Equation (3.8) shows that u and oj are linear coordinates on the hbers of D. Moreover, the vector 
helds Zi and Z 2 are invariant under the SE{2) action dehned in Section 2.3 and therefore they 


constitute a basis of sections of the reduced vector bundle D/SE{2). It follows that u and oj can be 
interpreted as linear coordinates on the hbers of the vector bundle D/SE{2) as advertised before. 


Equations (3.9) are of pure kinematic nature and are well known to the control community (see e.g. 
They dehne the evolution of the variables ai,... ,an in the reduced space and are consistent 


with (3.1). 


The evolution equation for u is of dynamical nature and can be easily obtained by noting that the 


nonholonomic constraints as written in (2.6) do not impose any restriction on the value of 9. Hence, 
the constraint reaction force written in the coordinates {x,y,9,9i, ... ,9n) has no component along 
the ^-direction, and the following dynamical equation holds 

d kdL\ dL ^ 


where L is given by (2.5). Explicitly we have 


(Jo + Ma^)'9 + Ma-^{y cos 9 — xsin^) = —Ma9{x cos 9 + ysm9). 


dt 


Using (3.9) we obtain 


OJ = 


Mauoj 
' Jo+ M a? 


(3.10) 


as in (3.1). 


The evolution equation for u is more difficult to obtain. As mentioned above, we follow the 
approach taken in [8]. This method to obtain the equations of motion of a nonholonomic system is 
outlined in the Appendix. 

The method requires us to compute the constrained Lagrangian Lc that is the restriction of L 
to H. It is the kinetic energy of the system when the nonholonomic constraints are satished. In 
view of the symmetries, its value can be expressed in terms of m, oj, ai,..., a„. To obtain an explicit 
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expression for Lc, start by noticing that (3.3), (3.7) and (3.9) imply 

fc-i 


X = u cos 9, y = usin 9, 9 = u, 9k = 


Next we prove the following. 


u sin ak 


n 


cos a,-, fc = 1,..., n. 


Proposition 3.2. Let j > 1. If the constraints (3.11) are satisfied, then we have 


x] + y] = cos^ Ofc. 


k=l 


Proof. By induction. The case j = 1 is a simple calculation using (2.4) and (3.11) and : 
reader. Assume that the result is valid for j — 1 > 1. Using (2.1) we write 


Xj = Xj-i + i9j sin 9j, ijj = — i9j cos 9j. 


Hence, 


2n2 


Xj + i/j = Xj_i + yj_^ + 2i9j{xj-i sin 9j - cos 9j) + £ 9^ 


Using (2.4) we write 


Xj-i = i; + £^sin 9k9k, i/j-i = y - i'^cos9k9k, 


k=l 


so that 


k=l 


1-1 


Xj_i sin 9j — ijj-i cos 9j = x sin 9j — y cos 9j + £ cos{9j — 9k)9k- 


k=l 


Now, in view of (3.11) and (3.3) we can write 

( ^ 

X sin 9j — y cos 9j = —u sin j 

\A:=1 




1-1 


'^cos(6'j - 9k)9, 


k — 


k=l 


+ n^cos 

k=l V^=A:-|-1 


^k-1 

ai I sm ak 1 cos as 


vS=l 


Using the identity (3.6) we conclude that 

Xj^i sin 9j — yj_i cos 9j = —£9j. 


Therefore, (3.13) becomes 


•2 I -2 -2 1-2 n 2 Q 2 

Xi+y, = x + 2/_1 -£ 9j. 


Using the induction hypothesis and (3.11) once more, this becomes 


1-1 


1-1 


• 2 I -2 2 TT 2 2-2 TT 2 

Xj + ijj = u M cos ak — u sm aj I I cos ak 


k=l 


k=l 


that is equivalent to (3.12). 


(3.11) 


(3.12) 

left to the 


(3.13) 
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It follows immediately from the above proposition, and from (3.11), that, if the nonholonomic 
constraints are satished, the kinetic energy %j of the trailer “Bj equals 


= d [JOj +m{x. +yA 


u 


0-1 

n 


\A;=1 


-K siy aj + m cos^ a,- 


for j = 1,..., n. For j = 0 we have 


^0 = ^ ( + Vc 


= - ((Jo + Ma^)u^ + Mu^) 


Therefore, adding up the contributions of all the cars in the convoy, we conclude that the con¬ 
strained Lagrangian is given by 


£jc — — (i?(ci)'U^ -|- (Jo -|- . 


(3.14) 


Next we prove the following. 


Lemma 3.3. The orthogonal projection of the commutator [y, Z 2 ] onto D with respect to the kinetic 
energy metric defined by the Lagrangian (2.5) is given by 


where 


PR{ay 


Q^Zi -l- 612.^2 

Q{^) ro 2 


p2 _ 
'-'12 ~ 


Ma 


Jo + Ma2 


Here Q{a) and R{a) are defined by (3.2). 


Proof. We give an indirect proof. In view of the discussion in the Appendix, the evolution equation 
for u can be obtained from the general formula (A.l) with the subindex 6 = 2 (for us = u, v'^ = u). 

by (P 

dL, 


Since Lc is independent of x, ?/, 9 and the vector held Z 2 is given by (3.4), we obtain 

d (dLr 


dt \ du 

where we have used C ®2 = — C 21 ) ^ 2. 


= eJsM 


du 


T CioH 




+ 


du dai ’ 


Using the expression (3.14) for Lc, the last equation becomes 


1 dR 


(Jo -|- Maf)u = + Ci 2 ('^o + Maf)uu ^ — — — {a)u^. 

2 do. I 


The above equation should simplify to (3.10) so we conclude that 


1 f)R 

e}2^(«) + 2^(«) = 0 , 


p2 _ 

'-'12 ~ 


Ma 


Jo + Ma2 


The proof is completed by noticing that 


dR 

dai 


[a = 


2Q(a) 
£2 • 


(3.15) 
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□ 


The equation for u can now be obtained from the general formula (A.l) with the subindex a = 1. 
Since Lc is independent of x, y, 9 and the vector field Zi is given by (|3.4l), we obtain 


^ / dLc \ _ dLc _ „2 9'^c 

dt[ du ) du du 




fe=i 


'dak 


that becomes 


On the other hand, 


d 1 ^ <5(«) 

j^(R{,o)u) = 


+ ^ ^ A} 


k=l 


OR 

da, 


-u 


d , T->f \ \ . nf \ ■ 

{R{a)u) = u } -—ak + Riam. 

oak 

k=l 


dt 


Using that 


we can combine (3.16) and (3.17) to give 


ai = oj + uAi, ak = uAk, /c = 2,..., r 

Q{a)uu 


R{a)u = 5^^/ 


,A :=1 


dR 

dak 


u 


dR 

dai 


UUJ 


Using ( |3.15| ) one shows that equation ( |3.19[ ) can be written as 

dR 


u = 


2R{a) 


v/c=l 


dak 


n" + 


Q(«) 

PR{a) 


P 


UOJ + 


T ]\d(ioj 


Ma 2 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


that completes the proof of Theorem m 


3.1. Energy conservation and the flow on the energy level surfaces. We note that, as it is 
usual with nonholonomic systems, the energy is preserved. In our case, this is the reduced kinetic 
energy given by the constrained Lagrangian (|3.14[). If we define 


£(«, u,u) = - [R{a)u^ + (Jo + Ma^)uP‘^ , 


(3.21) 


then a direct calculation that uses (|3.18|) and (|3.15|) shows that £ is preserved by the flow of equations 

(0. 

Let E > 0. It is natural to parametrize the level set 8, = E with the angles /3, oi,..., «„ where 
the angle f3 is uniquely determined by the conditions 


u = 


2E 

R{a] 


cos (3, 


u = 


2E 


Jo + Ma^ 


sin (3. 


(3.22) 


It follows that the energy level set £ = E is diffeomorphic to the {n + l)-torus To obtain an 

evolution equation for f3 we differentiate the above relation for u with respect to time to obtain 


U} = 


2E 


Jo + Ma2 


/3 cos (3. 
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Now, combining (|3.10|) with (|3.22|) and the above equation we obtain 

Ma 


2E 


Jo + Ma^ 


which simplifies to 


jS cos I3 = 


/S = 


2E 


Jo + Ma2 I ^R{a){Jo + Ma?) 


cos 13 sin (3 


Ma 


2E 


sin f3, 


(3.23) 


Jo + Ma'^ y R{a) 

assuming that cos/3 7 ^ 0. Proceeding in an analogous fashion, differentiating the relation for u in 
(3.22) with respect to time and using (3.20) we obtain (3.23) provided that sin/3 7 ^ 0. In conclusion, 
equation (3.23) holds for any value of /3. 

The rest of the equations for the flow restricted to the energy surface are obtained by combining 

-y/^sinai 

cos (3, 

(3.24) 


(3.22) with (3.1). We obtain 

sin /3 


(y.\ — 


ttfc = T 



2E 


Jo + Ma2 


i^yR{a) 


2E 

R{a) 


'k-2 


JJ^cosafc I (sinafc_i — coscLfc-i sinoffc) cos/3, 

\j=i 


k = 2,... ,n. 


We summarize the results of this subsection in the following. 


Theorem 3.4. The positive energy level sets of the reduced system ( 3.1[ ) are diffeomorphic to (u + 1)- 
tori that can be parametrized with the angular variables (/3,ai,... ,ari)- The restriction of the flow 
to the torus 8. = E > 0 is described by equations (3.23) and (3.24). 

4. Classification and linear stability of equilibria 

We study the equilibria of the reduced system restricted to a positive energy level set. Throughout 
this section we assume that the constant a > 0 . 

4.1. Classification of equilibria. 

Proposition 4.1. Let E > 0. There exist exactly 2”+^ equlibrium points in the energy level set 
8 = E of the reduced system (3.1). They are given by the conditions 

u = ±\l— -, (u = 0 , sinafc = 0 , fc = l,...,n. 


M + nm 


(4.1) 


Proof. We make use of the restricted equations (3.23) and (3.24). Imposing /3 = 0 in (3.23) implies 

sin/3 = 0. 


Under this condition, from (3.24) we see that we can only have di = 0 if sinoi = 0. Now assume 
that dfc = 0 and sin/3 = sinoi = • • • = sinafc_i = 0. From (3.24) it follows that sinofc = 0. This 
shows that the only equilibria of the system occur at the points where 

sin /3 = sin ai = ■ ■ ■ = sin = 0. (4.2) 


Now use (3.2) to show that the value of R{a) at these points is the total mass of the system 
M + nm. The proof is completed by using (3.22). □ 
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Assume that we are at an equilibrium configuration with energy E. The condition cj = 0 implies 
that the leading car moves along a straight line. It moves at the constant speed 


2E 

M+nm 


as indicated 


bygg, The motion is forward (in the direction from the midpoint of the wheel’s axis to the center 
of mass) if M > 0 or backwards if u < 0. 

On the other hand, the condition sinafc = 0 in (4.1) implies that the trailer 23^, is aligned with 
the {k — trailer "Bk-i- Denote by 

(Tfc = cos ttfc = ±1, fc = 1,..., n. (4.3) 

If CTfc = 1, then 23fc is ‘behind’ E>k-i- If cri = —1 then 23i ‘overlaps’ with 'Bq since we assume that the 
wheels of the leading car are located towards the rear of the vehicle. More generally, if cr^+i = — 1 
then 23fc_|_i ‘overlaps’ with Bk-i- See Figure 4.1 The situation resembles the equilibria of a chain of 
n coupled planar pendula. 


01=1 0^=1 011 + 1=1 
'Tj qj 'Xz q2 q2 

-^(1 -^1 -^k -^k+l 




Figure 4.1. Illustration of equilibrium states with ui = ±1 and with ak = 1 and (Tk+i = ±1. 

Therefore, the equilibria of the reduced system correspond to solutions where the convoy moves 
at constant speed along a straight line with all of the trailers aligned, with the possibility of overlaps 
between the cars. Of course the only physically attainable equilibria occur when cxi = ■ ■ ■ = = 1 so 

that there are no overlaps. There are two of such equilibria, corresponding to forward and backward 
motion of the convoy. We shall see that the former is asymptotically stable whereas the second one 
is asymptotically unstable. 

4.2. Stability of equilibria. We perform a linear stability analysis of the equilibria found in the 
previous subsection. We will consider the system restricted to the constant energy (n + l)-torus 
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E, = E, SO we work with equations (|3.23|) and (|3.24). To obtain the linearization of these equations 

(a) = 0, j = 


around an equilibrium, we shall use the relations 

OR 


R{a) = M + nm, 


dan 


that hold if a = (ai,..., a„) satishes the equilibrium conditions (4.1). 


Fix an equilibrium of equations (3.23) and (3.24) satisfying (4.2). Denote by 


(To = cos /3 = ±1. 

Forward motion of the convoy corresponds to (Tq = 1 and backward motion to (Tq = —1. 

A straightforward calculation shows that the constant (n + 1) x (n + 1) matrix that dehnes the 


linearization of (3.23) and (3.24) around the given equilibrium is 

/ 0 0 0 
-<^.<^0 0 0 


2E 


i \ M + nm 


0 

0 


(TiCTo 

0 


—(T2(Ti(To 
cr2(Ti(To 


0 

-(T3(T2(TiCro 


0 

0 

0 

0 




V 0 0 ■ ■ ■ ■ ■ ■ - n"=o / 

Since this matrix is lower diagonal, its eigenvalues are the diagonal components 


An — 


Ma 


2E 


Jo + Ma^ \ M + nm 


■Co, 


Afc — —T 


2E 


n 


a,- 


M + nm J-1 “ ’ 
i=o 


A; = 1,..., n. 


Therefore all of the equilibria are hyperbolic. Moreover, we immediately conclude the following about 
the nature of the equilibria. 

(i) If at least one of ak with A; = 1,..., n, is negative (there are overlaps between the trailers) 
then there are positive and negative eigenvalues and the equilibrium is a saddle point. 

If CTfc = 1 for all k = 1,... ,n (there are no overlaps) and (To > 0 (the convoy is moving 
forwards) then all of the eigenvalues are negative and the equilibrium is a stable node. 

If (Tfc = 1 for all k = 1,... ,n (there are no overlaps) and (To < 0 (the convoy is moving 
backwards) then all of the eigenvalues are positive and the equilibrium is an unstable node. 

An illustration of the numerical integration of the dynamics in the case n = 1 is given in Figure [T^ 
Here the constant energy surface is a two-torus. It is seen the the generic initial conditions approach 
the stable (respectively unstable) node as A —)• oo (respectively as A —)• —oo). Figure 4.2 also shows 


11 


111 


the trajectory of the leading car 'Bq on the plane for a generic initial condition. It asymptotically 
approaches steady motion along a straight line. The curve traced by Bq closely resembles the paths 
followed by the Chaplygin sleigh (see e.g. mm- 

5. The case a = 0 


If a = 0 the dynamics changes substantially. From (3.1) we see that u is constant throughout the 
motion. 
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(a) Phase portrait on a fundamental region of the torus (ai, /3). There (b) Trajectory of "Bq on the plane. Asymptotic 

are 4 equilibrium points (up to equivalence modulo 27r). A stable node behavior towards straight line motion, 

at (0,0), an unstable node at (0, tt) and two saddle points at (tt, 0) 
and (tt, tt). 


Figure 4.2. Phase portrait of the restriction of the reduced flow to an energy level two torus and 
generic trajectory of the leading car Sq in the case n = 1. 


If a; = 0, the classification of the equilibrium solutions of (3.1) coincides with the description given 


in Proposition 4.1, and the stability of the solution with 


(To — (Ji — ■■■ — (Jn — 1 


is analyzed in [7]. 

For the rest of the paper we consider the case where ca y 0. The classihcation of equilibria is more 
involved as the following proposition shows. 


Proposition 5.1. Suppose that a = 0 and that u = Uq ^ Q. A necessary and sufficient condition for 


the existence of equilibria of (3.1) with u = Uq is that 


nl 2 ^ 2 

ni (Uq < Uq. 


(6.1) 


Proof. Equations (3.18) imply that at such equilibria one must have uq ^ 0 and 


1^0 + uqAi — 0 , 


= 0, k = 2,... ,n. 


(5.2) 
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Using (3.5), the above equations can be written as 


iuJi 


sinoi = 


0 


cos sin 0^2 = 


Mo 

Mo 


n—1 

n 

k=l 


COS Q!fc Sin an = 


Mo 


One can inductively show that the solutions to the above equations satisfy 


Mq — 


£^uJq 


"" u^o-(k- "" ul-{k- l)eul ’ 

It follows that a necessary condition for the existence of equilibria is that 


k = 1,... ,n. 


ful 


Uq — {n — 1 )&Uq 


< 1 , 


which is equivalent to (5.1). That this condition is also sufficient is seen by noting that if (5.2) holds 


(and a = 0), the equation for u in (3.1) becomes 


. 1 fldR Q\ 

“ = s Uay + p I 


But the right hand side of this equation is zero by (3.15). 


□ 


The equations for x,y and 6 in (3.9) show that at an equilibrium solution with ^ 0 the car 


Bo moves along a circle of radius ^ at constant angular speed. Proposition (|5.1|) shows that the 


radius of this circle must be at least ^Jni. 

We do not attempt to study the stability properties of the system in this case. Instead, we treat 
the case of one trailer in detail. 


5.1. The case of one trailer. If a = 0 and n = 1 then, denoting ai = a, we have 

R{a) = M + m cos^ Oi ^ sin^ a, 


and the equations (3.1) become 


M = 


— J)u cos a sin a{luj — u sin a) 


i{{M + m cos^ a)£^ + J sim a) 


cu = 0, 


(6.3) 


Msina 


a = u — 


For physical reasons it is natural to assume 


J < m£^. 


(5.4) 
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Equations (5.3) are easily integrated using the conservation of energy. First notice that the level 


sets of the constants E and u are invariant circles parametrized by a 

u = 


12E — 


R{a) 


We £x a value of w = cuo > 0 and we study the behavior of the flow along the invariant circl^ 

(6.6) 


u = 


'2E — JquJq 


R{a) 


The evolution of a along the circle is given by 

R{a)uJo — \/2E — JoOOq sin a 


a = 


which leads to the quadrature 


R{a) da 


lujQ ^^/R{a) — \/2E — JquJq sin a 


= dt. 


(5.6) 


(6.7) 


Now notice that the inequality (5.4) implies 

J 


M + — < R{a) < M + m. 


Using (5.5) and the above inequality, we see that along the solutions of the system we have 


u‘ - Pull = - tVs 

R[a) 


< 


2E JqUJq 2 2 

•C Cc/n 


where 


M + ^ 

2(E- Ec 


Ec ■.= -{Jo + J + Mf) 


OJn. 


will depend on how E compares with Ec- 


The dynamics along the invariant circle 

Case 1. If ^JqUjI < E < E^.. 

It follows from Proposition 5.1 (or directly from (5.3)) that there are no equilibrium points of the 
system in this case. Hence, the dynamics along the invariant circle (5.5) is periodic. The energy 
dependent period T = T{E) is obtained using (5.7): 

I>2tt 

T = 


l\lR{a) da 


_ _=-. (5.8) 

Jo R{a) — ^y2E — JqUq sin a 

Using that E < Ec one can verify that the denominator does not vanish so this integral is convergent. 

Case 2. If E = Er. 


c 

'^The other cases, when either cjq or u or both are negative, are analogous. 
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In this case there is exactly one eqnilibrium point along the invariant circle (|5.5|) given by 

U = luJo, 


77 

cr = -. 


Hence, the invariant circle consists of a homoclinic connection and a critical point. 
Case 3. If > Ec- 
We have 


0 < 


< 


2E - Joul 2Ec - JquI J + MP' 


(5.9) 


The graph of the fnnction 


/(«) = 


sin^ a 


R{a) 


for 0 < a < 77 is shown in Figure 5.1 It is symmetrical with respect to a = 7r/2 where it achieves its 
maximum value of It attains every value between 0 and exactly two times. It follows 

from (5.9) that there exist exactly two values of a, that we denote by and a^‘^\ such that 

Pujo 


0 < < o < ^ ^^cl = 


2E — Jq^q 


2 ’ 


J = l,2. 



Figure 5.1. Plot of f{a). The horizontal line at height drawn under the assumption that 

E> Ec. 


A short calculation shows that the two points 


a = a 


U) 


sin(ahl) ’ 


J = l,2 


are the only equilibria of (5.3) contained in the invariant circle (5.5). 
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Given that sin(a('^)) > 0, j = 1,2, in a neighborhood of these points, we can write the evolution 
equation (5.6) for a as 

d = cuo - \j‘lE- f [a). 

Since / is increasing at and decreasing at we conclude that the equilibrium 


a = 


u = 


Uoi 


(5.10) 


sin(ab)) 

is asymptotically stable if j = 1 and asymptotically unstable if j = 2. A physical interpretation of 


these equilibria can be given with the aid of Figure 5.2 




Figure 5.2. Illustration of the unstable and stable equilibria if if > Ec- The configuration on the 
left is unstable and on the right is stable. The arrows indicate the motion of 130 (in agreement with 
our working assumption that u,uJo > 0). 


Hence, in this case, the invariant circle (5.5) consists of two heteroclinic orbits that connect the 
unstable critical point with the stable one. 

Figure 5.3 illustrates the qualitative dynamics on the invariant circle (5.5) in the three different 
energy regimes treated above. Figure |5^ illustrates the dynamics of (5.3) on the cylinder a; = cuo > 0. 



Figure 5.3. Qualitative behavior of the dynamics along the invariant circle (5.5) for the different 
energy regimes, E < aiv the left, E = Ec on the middle and E > Ec on the right. 
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a=0 


a 


Figure 5.4. Dynamics on the cylinder w = wq > 0. The vertical axis is u and the angular variable is 
a. It is seen that the motion takes place along the invariant circles. The critical points are indicated 
in black. The blue trajectories have E < Ec, the green ones E = Ec, and the red ones E > Ec- 


Our analysis shows that Ec is a critical value of the energy that separates two different qualitative 
behaviors. Subcritical energy values lead to periodic motion in the reduced space. On the other 
hand, supercritical energy values correspond to asymptotic behavior on the reduced space. A similar 
phenomenon is observed in the motion of a Chaplygin sleigh in a perfect fluid in the presence of 
circulation 


5.1.1. The motion on the plane. With the information given above, we can understand how the 2- 
body convoy moves in the plane. First note that in the absence of the trailer “Bi (i.e. if m = 0 and 
J = 0) then u = Uq, uj = loq for constants Uq and cuq. Hence, the motion of Bq on the plane for a 
generic initial condition is uniform circular motion on a circle of radius r = uq/\u:q\. 


Our analysis in the previous section shows that ii E > Ec in the limit as t —>■ ±oo the 2-body 
convoy on the plane approaches uniform circular motion. Continuing with the assumption that 
u^ujq > 0, from (5.10) we conclude that the radii of the limit circles is 



sin ad) sin ad) 


The value of sin ad) = sin ad) is decreasing and approaches 0 as the energy E ^ oo so the radius 
r —)■ cx) for large energies. Figure [575] shows a trajectory of the leading car obtained numerically. The 
trailer Bi locks itself at a hxed angle with respect to Bq as t —)■ ±cxd. The limit angles are ad) when 
t — )■ — CX) and ad) when f —)■ cx). 


On the other hand, if 0 < < Ec, the dynamics of a and u is periodic with period (5.8). After 

one period, the position of the leading car Bq suffers a rotation by an angle A6 = ujqT, followed by 
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(a) E < Eq. Generic 
quasiperiodic trajectory. 



(b) E < Ec- Periodic trajec¬ 
tory. 



(c) E > Ec. Asymptotic be¬ 
havior towards limit circles. 


Figure 5.5. Trajectory of the leading car iBg for different energy values. 


a translation by {Ax, Ay) with 

+ iAy = f u{t)e^‘^°^ dt = 


Ax 




'o ^JR{a{t)) 


dt, 


5.5 


where the dependence of a on t is determined by (5.7) and we have assumed that 6{0) = 0. 

Generically, the angle A6 is an irrational multiple of 27t and the motion of ‘Bq in the plane is 
quasiperiodic with its trajectory contained in an annulus or a circle. It is also possible to h ave 
periodic behavior if ^ G Q or unbounded trajectories if ^ G Z and Ax'^ + At/^ ^ 0. Figure 
shows a periodic and a quasiperiodic trajectory for Bq obtained numerically. 

6. Singular configurations 

The degree of nonholonomy is an important notion that arises in nonlinear control theory. It 
expresses the level of Lie-bracketing of the elements in the constraint distribution that is needed to 
span the tangent space at each conhguration. This concept comes up, for instance, when trying to 
quantify the complexity associated with steering the system from one point to another (see e.g. [TB] 1. 

When the number of trailers in our system is greater than or equal to two, this degree is not constant 


throughout the conhguration space. To hx ideas we treat the case n = 2 in detail. According to (3.4) 

/sin ai — cos ai sin 0^2 


a a sin ai a 

Zi = cos9— + sm9- - - —-— 

ux (jy t u(y.\ 


+ 


d 


d d 


Z2 — wx + 


\ i J da2 

form a basis of the constraint distribution D. Direct calculations show 

' cos «! -|- sin sin 02 


d9 dai ’ 


[Zi,Z 2 ] = sin 6 *^ - cos 9-^ + ^ 

OX uy t u(y.\ 


d 


[Zi,[Zi,[Zi,Z2]]] = 


cosui d 


1 d 


1 -|- cos 02 \ d 


2 -f- cos 02 


dao 


P dai 


COStti 


P 


A. 

da2 


dai ’ 

[Z2,[Zi,Z2]]=Zi, 

lZ2,lZi,lZi,Z2]]]=0. 


Let g G Q be a conhguration of the system with cosai 7 ^ 0. Then the vector helds Zi, Z 2 , [Zi, ^ 2 ], 
[Zi, [Zi, Z 2 ]] and [Zi, [Zi, [Zi, Z 2 ]]] form a basis of the tangent space TgQ. The element [Zi, [Zi, [Zi, Z 2 ] 
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in the basis is said to have length 4 since one needs to compute iterated brackets of four elements 
in the basis of 2) to generate it. It is clear that it is not possible to construct a basis for TqQ with 
iterated brackets of Zi and Z 2 and whose elements have length less than 4. We then say that the 
degree of nonholonomy at configurations with cos ai 7 ^ 0 is 4. 

On the other hand, at conhgurations q with coscti = 0, the vector held [Zi, [Zi, [Zi, Z 2 ]]] vanishes. 
One can complete Zi, Z 2 , [Zi, Z 2 ], [Zi, [Zi, Z 2 ]] to a basis of TgQ by adjoining the vector held 


[[Zi,Z2],|Z„|Zi,Z2]]] 


sinai d , /2 + cosa2\ d 

P dai ~ V 


that has length 5. Hence, the degree of nonholonomy at conhgurations with cosoi 7 ^ 0 is 5. 

The latter conhgurations are called singular and correspond to having ‘Bq jackknifed, that is, Bq 
and Bi are perpendicular. It is intuitively clear that maneuvering the system at this conhguration 
is a more difficult task. The classihcation of singularities for the n-trailer vehicle, and the degree of 
nonholonomy at each of them, is given in [U] for arbitrary n. These correspond to diherent jackknihng 
possibilities for the bodies in the convoy. A natural question is to understand what are the ehects of 
these singular conhgurations on the dynamics, if any. 

Another example of a nonholonomic system exhibiting singular conhgurations is an articulated 
arm. In recent years there have been diherent ehorts to classify the singularities of the associated 
constraint distribution [mil]. 

To our knowledge, the ehect of this kind of singularities on the motion of nonholonomic systems 
is unexplored. We hope to report on this issue in a future note. 


Appendix. 


The derivation of the evolution equation for m in (3.1) relies on the method given in [8] to obtain the 
equations of motion of a mechanical nonholonomic system. This reference includes a more detailed 
description of the geometry and considers more general cases than what we need. Here we only outline 
the main steps to obtain (a simple version of) their equations (3.7) and (3.8). Our presentation is 
done without proof. 


Consider a nonholonomic system on a conhguration manifold Q of dimension N with Lagrangian 
£ : TQ —)■ M of mechanical type and constraint distribution B) of constant rank k < N that is bracket 
generating. The condition that £ is of mechanical type means that it is the sum of kinetic minus 
potential energy, and that the kinetic energy dehnes a Riemannian metric 9 on Q. 

Associated to the metric 9 there is a decomposition TQ = B) ® R-*-, where is the 9-orthogonal 
complement of R, and a projection CP : TQ — )• R. 

The idea is to write down the equations of motion that are consistent with the Lagrange d’Alembert 
principle using quasi-velocities that are adapted to the distribution R. Denote by , local 

coordinates on an open set of Q and by Zi, ... ,Zk a basis of sections of R in such open set. That 
is, they are linearly independent vector helds that lie on R. 

Dehne the scalar functions pi and on Q through the relation^ 

= ‘^{.[Zb^Zd]) = Ql^Ze, 6 ,d,e = 1,.../c, i = l,...,N. 


6 


Here and in what follows we use the convention of sum over repeated indices 
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Let q G Q. Any tangent vector n G 2)^ C TqQ can be written as 

V = v'^Zbiq) 

for certain scalars (the quasi-velocities). Hence, the value of the restriction of the Lagrangian to 
H, that we denote as can be expressed in terms of the variables q^,... ,q^,v^,... ,v^. 

Equations (3.7) and (3.8) in [8] state that the equations of motion for the nonholonomic system can 
be written as 


q" = P\v\ 
d / dLc\ 
dt \ dv^) 


rif 

Pi 


dv^ 


dq^ ’ 


6 = 1,..., fc. 


(A.l) 


These equations avoid dealing with Lagrange multipliers. The effect of the constraint forces is 
encoded in the effect of the projector T on the dehnition of the structure coefficients 
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